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Abstract 

The AMIDAS website has been estabhshed as an onUne interactive tool for running simu- 
lations and analyzing data in direct Dark Matter detection experiments. At the first phase 
of the website building, only some commonly used WIMP velocity distribution functions 
and elastic nuclear form factors have been involved in the AMIDAS code. In order to let 
the options for velocity distribution as well as for nuclear form factors be more flexible, 
we have extended the AMIDAS code to be able to include user-uploaded files with their own 
functions. In this article, I describe the preparation of files of user-defined functions onto 
the AMIDAS website. Some examples will also be given. 



1 Introduction 



In the last few years we developed new methods for analyzing data, i.e., measured recoil energies, 
from (future) direct Dark Matter detection experiments as model-independently as possible 
[1, 2, 3, 4]. These methods will help us to understand the nature of WIMP (Weakly Interacting 
Massive Particle x) Dark Matter, to identify them among new particles produced hopefully in 
the near future at colliders, as well as to reconstruct the (sub) structure of our Galactic halo. 
Following the development of these model-independent data analysis procedures, we combined 
the programs for simulations to a compact system: AMIDAS (A Model-Independent Data Analysis 
System). For users' convenience and under the collaboration with the ILIAS Project [5], an online 
system has also been established at the same time [6, 7]. 

For the first version of the AMIDAS code and website, the options for target nuclei, for the 
velocity distribution function of halo WlMPs, as well as for the elastic nuclear form factors 
for spin-independent (SI) and spin-dependent (SD) WIMP-nucleus interactions are fixed and 
only some commonly used forms have been involved in the AMIDAS code [7]. Users can not 
choose different detector materials nor use different velocity distribution/form factors for their 
simulations and/or data analyses. In order to let the options for the velocity distribution as well 
as for the nuclear form factors be more flexible, we have extended the AMIDAS code to be able to 
include user-uploaded with their own functions. Note that, since the AMIDAS code has been 
written in the C programming language, all user-defined functions for uploading must also be 
given using the syntax of C. 

The remainder of this article is organized as follows. In Sec. 2 1 will talk about setting users' 
own target nuclei. In Sees. 3 and 4 the preparation of files defining the velocity distribution func- 
tion and the nuclear form factors will be described, respectively. I will conclude in Sec. 5. Some 
intrinsically defined constants and functions in the AMIDAS code will be given in an appendix. 



2 Target nuclei 

Users can give as usual the element symbol and the atomic mass number A of their chosen 
target nuclei on the website directly. For the determination of ratios of different WlMP-nucleon 
couplings/cross sections, the total spin of the target nucleus J and the expectation values of the 
proton and neutron group spins (>S'(p^n)) for the chosen targets with spin sensitivities are also 
required. 



3 Velocity distribution of halo WIMPs 



For defining the one-dimensional WIMP velocity distribution, users should in practice give the 
integral over the velocity distribution function: 



/■fmax 



fiiv,t) 



dv 



(1) 



with the name of "Intf lv_v_user". Note here that the lower limit of the integral, the minimal 

incoming velocity of incident WIMPs that can deposit the energy Q in the detector, fmin, must 
be expressed as a function of the energy Q through f min = OL\fQ] the upper limit of the integral, 
^^max, is often set as the escape velocity fesc or even cxo, since the WIMP fiux on the Earth is 



usually assumed to be negligible at velocities v > Vesc ^ 600 km/s. Meanwhile, since 



a = 



mN 



r,N 



with the reduced mass 

mr,N = \ 



(2) 



(3) 



is a function of the WIMP mass and the mass of target nucleus mN, which has been defined as 
a function of the atomic mass number A, the expression of the integral over the WIMP velocity 
distribution should also be a function of m^ and A. Finally, the velocity distribution function 
and thus the expression of the integral should generally be a function of time t. 

As examples, the integral over the simple Maxwellian velocity distribution function [8, 1] can 
be defined as 

Example 



double Intflv_v_user (double mchi, int A, double QQ, double tt) 
{ 

return 

(2.0 / sqrt(M_PI) / (v_0 * v_U)) * 

exp (-alpha (mchi , A) * alpha (mchi , A) * QQ / 

( (v_0 * v_U) * (v_0 * v_U) ) ); 

> 



and the integral over the shifted Maxwellian velocity distribution [8, 1] can be given by 
Example 



double Intflv_v_user (double mchi, int A, 
{ 


double QQ, 


double tt) 




return 








(1.0 / 2.0 / (v_e(tt) * v_U)) * 








( erf( (alpha (mchi. A) * sqrt(QQ) + 


(v_e(tt) * 


v_U)) / (v_0 


* v_U) ) 


- erf( (alpha (mchi, A) * sqrt(QQ) - 

> 


(v_e(tt) * 


v_U)) / (v_0 


* v_U) ) ) ; 



Here v_e(tt), standing for the time-dependent Eaith' s velocity in the Galactic frame Ve{t), has 
been defined in the AMIDAS code by 



Definition 



double v_e (double tt) 
i 

return v_0 * ( 1.05 + 0.07 * cos (omega * (tt - t_p)) ); 

} 



Note that, firstly, v_U, the function alpha (mchi, A), where mchi and A stand for the WIMP 

mass m^ and the atomic mass number of the target nucleus. A, and the constant omega are 
defined intrinsically in the AMIDAS code. Secondly, v_0 and t_p, standing for the Sun's orbital 
velocity in the Galactic frame, Vq, and the date on which the velocity of the Earth relative 
to the WIMP halo is maximal, tp, are two input parameters which users can set later on the 
website [7]. Thirdly, QQ and tt stand for the energy variable Q and time variable t. Finally, 



for a time-independent velocity distribution or a numerical expression with a fixed experimental 
running time t — texpt, "double tt" should still be declared as one of the /our variables of the 
function Intf lv_v_user. 

On the other hand, for comparing results of the reconstructed velocity distribution with 
the (input) theoretical one in the output plots, users need an extra file for giving the velocity 
distribution function itself (not the integral over it now) for drawing the (input) theoretical 
velocity distribution function fi{v). Since the gnuplot software [9] has been used for drawing 
output plots, the velocity distribution function given in this file must be written using the syntax 
of gnuplot with the name of "f lv_user(x)". As examples, the simple Maxwellian velocity 
distribution can be given as 

Example 

M_PI = 3.1416 

v_0 = 220.0 

flv_user(x) \ 

= (4.0 / sqrt(M_PI)) * \ 

((x * x) / (v_0 * v_0 * v_0)) * \ 

exp(-(x * x) / (v_0 * v_0)) \ 

and for the shifted Maxwellian velocity distribution: 
Example 

M_PI = 3.1416 

omega = 2.0 * M_PI / 365.0 

t_p = 152.5 
t_expt = 243.75 

v_0 = 220.0 

v_e = v_0 * ( 1.05 + 0.07 * cos (omega * (t_expt - t_p)) ) 

flv_user(x) \ 

= 1.0 / sqrt(M_PI) * (x / v_e / v_0) * \ 

( exp(-(x - v_e) * (x - v_e) / (v_0 * v_0)) \ 

- exp(-(x + v_e) * (x + v_e) / (v_0 * v_0)) ) \ 

Note here that the "\" (backslash) must be used (also at the end of the last line!) in order to 
insert the definition given in this file into the other intrinsic commands correctly. 

Two sample files, one is for the AMIDAS code and the other one is for the gnuplot software, 
can be downloaded from the AMIDAS website. 

4 Nuclear form factors 

For defining users' own nuclear form factors, not only the definition of the squared form factor 
F'^{Q) but also its derivative with respect to the energy Q: 



dF^ 



iQ) 



2F{Q) 



dF{Q) 
dQ 



(4) 



dQ 



(not F{Q) itself nor dF{Q) / dQ) must be given in one file together with the names of "FQ_SI_user 
(FQ_SD_user)" and "dFC)dQ_SI_user (dFQdQ_SD_user)" . The expression of the squared form 
factor can be either a general analytic form for different nuclei or a specified form for the cho- 
sen nucleus. As examples, for the SI WIMP-nucleus cross section, the exponential form factor 
F^^{Q) [10, 11, 8] can be given as 

Example 

double FQ_SI_user(int A, double QQ) 
i 

return exp(-QQ / Q_0(A)); 

> 

double dFQdQ_SI_user(int A, double QQ) 
{ 

return -(1.0 / Q_0(A)) * FQ_SI_user (A, QQ) ; 

} 



and the Woods-Saxon form factor F^q{Q) [12, 8]: 



Example 



double FQ_SI_user(int A, double QQ) 








{ 








if (QQ == 0.0) 








return 1.0; 








else 








{ 








return 








( 3.0 * sphBesselj (1, qq(A, QQ) * R_1(A)) / 


(qq(A, 


QQ) 


* R_1(A)) ) * 


( 3.0 * sphBesselj (1, qq(A, QQ) * R_1(A)) / 


(qq(A, 


QQ) 


* R_1(A)) ) * 


exp(-(qq(A, QQ) * ss) * (qq(A, QQ) * ss)); 








} 








> 








double dFQdQ_SI_user(int A, double QQ) 








{ 








if (QQ == 0.0) 








return -0.4 * (R_A(A) * R_A(A)) * in_N(A); 








else 








{ 








return 








( dsphBesselj (1, qq(A, QQ) * R_1(A)) * R_1(A) / 






sphBesselj (1, qq(A, QQ) * R_1(A)) 








- 1.0 / qq(A, QQ) 








- (qq(A, QQ) * ss * ss) ) * 








FQ_SI_user(A, QQ) * ((2.0 * m_N(A)) / qq(A, 


QQ)); 






} 








> 









Here the spherical Bessel functions 




(5) 



their derivatives 




^+1/2 (a;) + 4+1/2(2^) > 



(6) 



as well as the (derivatives of the) half-integer Bessel functions J„+i/2(a;), for n = 0, ±1, ih2, ■ ■ ■, 
are given in an intrinsic package of the AMIDAS code. The functions Q_0 (A) , qq(A , QQ) , R_l (A) , 
R_A(A), m_N(A), and the constant ss are also defined intrinsically in the AMIDAS code (see the 
appendix) . 

As a more complicated example, the thin-shell form factor i^Ts(Q) used sometimes for the 
SD WIMP-nucleus cross section [13, 14] can be defined as 

Example 

double FQ_SD_user(iiit A, double QQ) 
{ 

if (QQ == 0.0) 
return 1.0; 

else 

if (QQ <= QQ_SD_min(A) | | 
QQ >= QQ_SD_max(A) ) 

■C 

return 

sphBesselj (0, qq(A, QQ) * R_1(A)) * sphBesselj (0, qq(A, QQ) * R_1(A)); 

} 

else 

return FQ_TS_const; 

> 

double dFQdQ_SD_user(int A, double QQ) 
{ 

if (QQ == 0.0) 

return -(1.0 / 3.0) * (R_1(A) * R_1(A)) * (2.0 * m_N(A)); 
else 

if (QQ <= QQ_SD_min(A) | | 
QQ >= QQ_SD_max(A) ) 

{ 



dsphBesselj (0, qq(A, QQ) * R_1(A)) * R_1(A) * 

sphBesselj (0, qq(A, QQ) * R_1(A)) * ((2.0 * m_M(A)) / qq(A, QQ)); 

} 

else 

return 0.0; 

} 



return 



The functions QQ_SD_min(A), QQ_SD_max(A), and the constant FQ_TS_const are defined in the 
AMIDAS code (see the appendix). 

5 Summary 

In this article, I described the preparation of files giving user-defined WIMP velocity distribution 
function and/or elastic nuclear form factor (s) which can be uploaded onto the AMIDAS website 
for more flexible simulations or data analyses. This improvement allows theorists to simulate 
with their own/favorite models and compare them with (future) experimental results, as well 
as gives experimentalists flexible choices for more suitable form factor (s) for their own detector 
materials. 

In summary, up to now all basic functions of the AMIDAS code and website have been well 
estabhshed. Hopefully this new tool can help our colleagues to detect/discover WIMP Dark 
Matter, to understand the nature of Dark Matter particles and the (sub) structure of the Galactic 
halo in the future. 
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A Relevant constants and functions defined in the AMIDAS 
code 

Here I list some relevant constants and functions defined intrinsically in the AMIDAS code, which 
could be needed by users for defining their own functions. 



A.l Defined constants 

The constants defined in the AMIDAS code are in natural units as following: 
Definitions 



double 


c = 1.0; 




double 


v_U = c / (2.998 * 


le5); 


double 


m_U = le6 / (c * c; 


; 


double 


fm_U = c / (0.1973 


* le6); 


double 


m_p = 0.938 * m_U; 




double 


omega = 2.0 * M_PI 


/ 365.0; 



A. 2 Defined functions 

The functions used for defining the commonly used WIMP velocity distribution functions and 
nuclear form factors in the AMIDAS code are given here. 



A. 2.1 The reduced mass mr,N(»Tij^, A) and cx{my^, A) 

All masses in the AMIDAS code are in the unit of GeV/c^. Firstly, the mass of the target nucleus 
is defined as a function of the atomic mass number A by 



Definition 



double m. 


.N(int 


A) 


{ 






return 


m_p * 


A * 0.99; 


> 







Here the mass difference between a proton and a neutron has been neglected. And the reduced 
mass between the WIMP mass and "something" is given by 

Definition 

double mchi_r (double mchi, double mx) 
{ 

return mchi * mx / (mchi + mx) ; 

> 



Therefore, the reduced mass of the WIMP mass with the mass of the target nucleus, mr,N in 
Eq.(3), and with the proton mass are defined as 



Definitions 



double mchi_rN (double mchi. 


int A) 


{ 




return mchi_r(mchi * m_U, 


m_N(A)); 


> 




double mchi_rp (double mchi) 




{ 




return mchi_r(mchi * m_U, 


m_p) ; 


> 





Finally, a defined in Eq.(2) has been given as a function of and A: 

Definition 

double alpha(double mchi, int A) 
{ 

return sqrt(m_M(A) / 2.0) / mchi_rN(mchi , A); 

> 



A. 2. 2 For nuclecir form factors 

For the exponential form factor Fex{Q), the nuclear radius Rq and the nuclear coherence energy 
Qo have been defined as functions of the atomic mass number A: 

Definitions 

double R_0(int A) 
{ 

return (0.3 + 0.91 * cbrt(m_N(A) / m_U)) * fm_U; 

> 



double Q_0(int A) 
{ 

return 1.5 / (in_N(A) * R_0(A) * R_0(A)); 


For the Woods-Saxon form factor -Fws(Q)) the nuclear skin thickness s, the effective nuclear 

radius Ra, and the radius i?i are defined as 

Definitions 

double ss = fiii_U; 

double R_A(int A) 

{ 

return 1.2 * cbrt(A) * fm_U; 

> 

double R_l(int A) 
{ 

return sqrt(R_A(A) * R_A(A) - 5.0 * ss * ss) ; 



Meanwhile, the transferred 3-momentum q = \/2rn^Q has been give as a function of the atomic 
mass number A and recoil energy Q: 

Definition 

double qq(int A, double QQ) 

{ 

return sqrt(2.0 * m_N(A) * QQ) ; 
_]_ 

For the thin-shell form factor Fts(Q), the dimensionless constants giving the lower and upper 
energy bounds, between which the form factor is a constant, and this constant (~ 0.047) are 
given as 

Definitions 

double qqR_l_min = 2.55; 
double qqR_l_max = 4.50; 

double FQ_TS_const ; 

FQ_TS_const = sphBesselj (0, qqR_l_min) * sphBesselj (0, qqR_l_min) ; 



Then the lower and upper energy bounds can be estimated by 
Definitions 



double QQ_SD_min(int A) 














i 














return (qqR_l_min * qqR_ 


l_min) 


/ 


(2.0 


* m_M(A) 


* R_1(A) 


* R_1(A)); 


> 














double QQ_SD_max(int A) 














i 














return (qqR_l_max * qqR_ 


l_max) 


/ 


(2.0 


* m_M(A) 


* R_1(A) 


* R_1(A)); 


} 















A. 3 Declared variables 

Here I list the short used names for variables and constants in the AMIDAS code. Note that these 
names should be avoided to use in the user-uploaded code(s), otherwise the AMIDAS program 
might not work correctly. 

Definitions 

AX, AX_str, Ax, calN, fp_U, G_F, hn, InX, JJ, k_l, ka, kn, lambda. In, 
mchi_assiimed, mchi_tmp, mm, m_N_TX, nexpt , nmchi , NmX, nn, iin_max, nn_win, 
nplot, nranap, nTX, pb_U, pm, Qmax, Qmid, Qmin, Qwin, Qvin, Qtp, rho, rho_0, 
rho_U, RJ, RJX, rlh, Rn, Rii_sh, RnX, Rsigma, R_TX, run, runi, runj , runk, 
r_win, sh, Sn, Snp, sol, Sp, Spn, std_rho_0, std_t_p, std_v_0, s_U, t_end, 
t_expt, tmax, tmid, tmin, t_start, TX, v_e_annual, v_esc 
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